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We present branching fraction measurements of twelve B meson decays of the form B — » 
£)(,*) D^jy The results are based on X(4S) decays in BB pairs. One of the B mesons is fully recon- 
structed and the other decays to two charm mesons, of which one is reconstructed, and the mass and 
momentum of the other is inferred by kinematics. Combining these results with previous exclusive 
branching fraction measurements, we determine B(D^ — > 4>n~) = (4.62 ± 0.36 s tat. ± 0.51 sys t.)%, 
B(D s j(2460)~ -> D* a ~ir°) = (56 ± 13 stat . ± 9 syst .)% and B(D s j(2A60y -> £>J 7 ) = (16 ± 4 sta t. ± 

3syst. )%• 

PACS numbers: 13.25.Hw, 13.25.Ft 



In this paper we present the study of charged and 
neutral B mesons decaying to two charm mesons, i.e. 
B — > D meas Dx U- -Dmeas represents a fully reconstructed 
£)(*)+>Q or £)(*) meson, and the mass and momentum 
of the Dx are inferred from the kinematics of the two- 
body B decay. This study allows measurements of B 
branching fractions without any assumption on the de- 
cays of the Dx ■ Measurements of these two-body branch- 
ing fractions can provide tests of the factorization of 
the decay amplitudes in the high momentum trans- 
fer regime 3]. From two separate classes of events with 

(*)— (*)— 
D 'mens — Ds and with Dx = D s we measure the 

branching fraction of Dj — > 4>7r~ , which has important 
implications for a wide range of D s and B physics. Fur- 
thermore, we select final states with Dx = D s j (2460)" 
and combine with the BABAR measurements of B(B — > 
£>(*)+>Q.D sJ (2460)-) x S(£> s j(2460)- -> Z?*"^ ) and 
B(B -> D(*)+'°L> sJ (2460)") x B(D sJ (2A60)- -> 
D~j) 4], thus extracting for the first time the absolute 
branching fractions of this recently observed state || . 

This analysis uses T(4S) — > BB events in which either 
a B + or a B° meson decays into a fully reconstructed 
hadronic final state (B r0 co)- The measurements are based 
on an integrated luminosity of 210.5 fb" 1 recorded at 
the T(4S) resonance with the BABAR detector at the 
PEP-II asymmetric-energy e + e~ collider operating near 
the T(4S) resonance. An additional 21.7 fb _1 recorded 
40 MeV below the resonance (off-resonance) are used to 
evaluate backgrounds. The BABAR detector is described 
in detail elsewhere 0. Charged-particle trajectories are 
measured by a vertex tracker with 5 double-sided layers 
and a 40-layer drift chamber, both operating in a 1.5-T 
magnetic field of a superconducting solenoid. Charged- 
particle identification is provided by the specific energy 
loss (dE/dx) in the tracking devices and by an internally 
reflecting ring-imaging Chcrcnkov detector. Photons are 
detected by a CsI(Tl) electromagnetic calorimeter. Wc 
use Monte Carlo simulations (MC) of the BABAR detec- 
tor based on GEANT4 [7J to optimize selection criteria 
and determine selection efficiencies. 

To reconstruct a large sample of B mesons, the 



hadronic decays B rcco — * DY + ,D*Y + are selected. 
Here, the system Y + consists of hadrons with a to- 
tal charge of +1, composed of ni7r ± n-iK^ n^K® 77,471"°, 
where tlx + 712 < 5, n% < 2, and 71,4 < 2. We re- 
construct D*~ -> D°ir-; D*° -> D°ir ,D T , D~ -> 

K+1T-TT-, K+1T-TT-1T , K®TT~ , K°TT~ 7T° , K°TT~ 7T" 7T+ J 

D° -> K+7T-, K+n-ir , K + ir-K-iT + , K° s it + it~- and 
Kg — > 7r + 7r". The kinematic consistency of B rcC o can- 
didates is checked with two variables, the beam energy- 
substituted mass toes = y/s/i—pg and the energy dif- 
ference AE = Eb —y/s/2. Here y/s is the total energy in 
the T(AS) center-of-mass (CM) frame, and %>b and Eb 
denote the momentum and energy of the B reco candidate 
in the same frame. The resolution on AE is measured to 
be (Tab = 10 — 35 MeV, depending on the decay mode, 
and we require \AE\ < 3cja.e- 

For each reconstructed B decay mode, the purity V is 
estimated as the ratio of the number of signal events with 
toes > 5.27GeV/c to the total number of events in the 
same range, and is evaluated on data. We only use modes 
for which V exceeds a decay-mode dependent threshold 
in the range of 9% to 24%. In events with more than one 
Brcco we select the decay mode with the highest purity. 
On average, we reconstruct one signal B rcco candidate in 
0.3% (0.5%) of the B°B° (B + B~) events. 

The selected sample of B rcC o is used as normalization 
for the determination of the branching fractions. It is 
contaminated by e + e~ — > qq (q = u,d,s,c) events and 
by other T(45) -> B°B° or B+B~ decays, in which the 
B roC o is mistakenly reconstructed from particles coming 
from both B mesons in the event. To significantly reduce 
the e + e~ — > qq background we require the angle 
defined in the CM frame, between the thrust axis |8| of 
the B loco and the thrust axis of all charged and neutral 
particles in the event excluding the ones that form the 
-B re co j to satisfy the requirement | cos8^ B \ < 0.7. 

On this signal-enriched sample (Fig. Q), the contribu- 
tions from the background are estimated as the sum of 
three components: the e + e~ — * qq, the B°B°, and the 
B + B _ events. The shapes of these background distribu- 
tions arc taken from MC simulation. The normalization 
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FIG. 1: Distributions in tubs for the -B rGC o sample. The back- 
ground contributions, determined as described in the text, are 
overlaid. 



of the e + e~ — > qq background is taken from off-resonance 
data, scaled by the luminosity. The normalization of the 
B a B°, B + B~ components are instead obtained by means 
of a x" 2 fit to the toes distribution in the sideband region 
(5.21GeV/c 2 < toes < 5.26GeV/c 2 ). The background 
contamination in the signal region (toes > 5.27GeV/c 2 ) 
is extrapolated and subtracted from the data to esti- 
mate the signal yield. After correcting for the | cos8'^ B \ 
cut efficiency estimated in the MC, the size of the to- 
tal sample of fully reconstructed B decays is N B o — 
(2.90±0.01 st at.)xl0 5 andTV B + co = (4.63±0.01 sta t.) x°10 5 . 

From the charged tracks and the neutral clusters that 
do not belong to the B reco we reconstruct the charmed 
mesons (D mea s) m the modes D° — > K~ir + , K~ir + ir°, 
K-TT+TT-TT+; D+ -> K-TT+TT+, D+ -> K°ir+ ; and 
D~ -> (f>Tr- (<f> -> K+K~), K° S K- {K° s -> vr+vr-), 
and K*"K- (K*° -> K+n~). We select and K*° 
candidates with a reconstructed mass within 15MeV/c 2 
and 70MeV/c 2 from their nominal values 0, respec- 
tively. The D* candidates are reconstructed in the decay 
modes D*+ -> D tt+, D+ir°, D*° -> D°ir°, L> 7 , and 
D* s ~ — ► -D77. We require the reconstructed masses of 
the D°, D + , and D ~ candidates and the differences Am 
between the masses of the D* and D candidates to be 
within 1.5 — 3 times its measured resolution from their 
nominal values [||, depending on the background level. 

We apply further selection criteria to enhance the 
signal contributions in the sample. For D^ + D x and 

Dg D x we consider neutral i? rcC o candidates while for 
£)(*)0£)- anc j £>(*) £>o_ we re q U j re positive charged -B reC o 
candidates. We suppress background from B — > D^*Hv, 
while keeping events with a semileptonic Dx decay, by 
rejecting any event with a remaining identified lepton 
with the appropriate charge and a momentum in the B 
rest frame (p*) greater than lGeV/c. In order to mini- 
mize the contamination of the modes with a D* to the 
modes with a D meson, we assign the events consistent 
with both the hypotheses (B — > DD X and B — > D*D X ) 
to the D* sample. 



The invariant mass of Dx {mx) is derived from the 
missing four- momentum p x = Pt{as) — Pb icco - Po meas , 
where all momenta are measured in the laboratory frame. 
The mx resolution is improved by a global Y(4S) kine- 
matic fit 01 that includes beam position and energy in- 
formation and constrains the masses and decay vertices 
of the -Dmeas- The x 2 of this fit is used to reduce the 
combinatorial background. We remove reconstructed D 
mesons with x 2 probability smaller than 0.1%. 

Of the selected events, 3 — 6% (9 — 30%) contain multi- 
ple D( s ) (Dfs) candidates. We retain those in the D mcstS 
decay mode with the lowest combinatorial background. If 
there are multiple candidates with the same decay mode, 
we select the one with the lowest value of |tod — mpDcl 
and K^-mpflG) 2 /^ + {Am-Am PDG ) 2 /a 2 Am 
for Z?( s ) and Dt^ respectively, where to is the recon- 
structed mass of the -D mea s candidate and the subscript 
PDG indicates nominal values 0. 

Finally, we consider only candidates in the range 
1.65GeV/c 2 < m x < 2.71 GeV/c 2 for the D^+/°D X 
modes and 1.68GeV/c 2 < m x < 2.31 GeV/c 2 for 
D y s '~D X - These ranges were chosen to minimize the to- 
tal uncertainty introduced by the background shape and 
normalization. 

The yield of each decay mode is extracted from the mx 
distribution by a binned x 2 fit of a sum of n s - lg signal con- 
tributions (N S1S ) and the total background contribution 
(7V bkg ), which is a sum of the combinatorial background, 

(*) (*) 
other B — > D\ID X decays, and D)t — crossfeed, 

to the experimental data. The signal and background 

distributions are histograms taken from MC simulation. 

For D°Dx we also weight the background shape with a 

second order polynomial function whose parameters are 

fitted on data. In the case 

of D^+/°D X modes we con- 
sider three signal components: D^ + /°D~, Z)(*)+/ Z)*-, 
and L>W+/°.D s ,/(2460)-, while in the case of D ( * ] ~D X 
modes we consider two signal components: dI*^ Z?+/° 
and D^~D*+/°. The x 2 is defined as: 

^ir r \ \ " I AT lcas ~ /^(Cj,Cbk g ) \ 
i \^/5N™ cas2 + 6N? IC J 

where 7V™ oas is the number of observed events in bin i, 

Hi corresponds to fa = Y^j=i, n , ig C i N if + Cbkg^ ks , 
the index j denotes the signal component, and <5A^ meas 
and SN^ IC are the statistical uncertainties for data and 
MC samples, respectively. The relative normalizations of 
each component (Cj and Cbkg) are allowed to vary in the 
fit. The measured mx distributions and the results of 
the fits are shown in Fig. [2J 

The branching fractions are extracted as B(f) = 
Nfi t /(sNB rcoo ), where Nfn is the number of signal events 
obtained from the fit to the mx distribution for a given 
mode and e, which includes the intermediate branching 
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FIG. 2: Distributions of m x . Fitted B -> D M+ '° and B -> D ( * )+ ' L> s j(2460)~ signal contributions and background 
components, determined as described in the text, are overlaid to the data points. 



fractions of -D moa s and its decay products, is the selection 
efficiency estimated using MC simulation. 

The dominant systematic uncertainties originate from 
the lack of knowledge of the correct shapes used in the 
mx fit, and from the determination of efficiencies (be- 
cause of the limited MC statistics). These uncertain- 
ties range from 5.6% to 25%, depending on the mode. 
The systematic uncertainties due to the determination 
of Nb ibco and to the differences between data and MC in 
the composition of the reconstructed -B reC o modes range 
between 3.7% and 6.7% for B°, and between 3.5% and 
9.0% for B + depending on the mode under study. Other 
uncertainties come from track reconstruction efficiency 
(1.4% per track and 2.2% per soft pion), 7 and 7r° effi- 
ciencies (3.0% per 7T and 1.8% per 7), and kaon identifi- 
cation (2% per kaon). The uncertainties due to branching 
fraction measurements for exclusive D^*} decays con- 
tribute between 3.0% and 7.4%, depending on the mode. 
We check the uncertainties introduced by the \ 2 cu t of 
the kinematic fit by comparing data and MC control sam- 
ples for B — > D^*'lv obtained with all previously men- 
tioned cuts except for the p* > 1 GeV/ c criterion applied. 
The statistical uncertainty of this comparison is used as 
the systematic uncertainty (between 0.5% and 2.3%). 

We combine the sixteen measurements of B — > 
£)(*)£)(*) to obtain the eight branching fractions for 

these modes and B(D~ — * 4>ir~) in a \ 2 fit. In this com- 



bination the ratios B(D~ -> K*°K-)/B{Dj -> (f>n~) 
and B(D- -> K°K~)/B(Dj -► <fm~), included in the 
efficiency calculation when -D m cas = D*^ , are fixed 
while B(D~ — > 4>tt~) is a free parameter. The MC model 
used to generate the D~ — > K + K~n~ decays does not 
include any interference among the different final states 
(07r~, K*°K~ , /o(980)7r _ , ...). Correlated and uncorre- 
lated uncertainties are properly taken into account in the 
covariance matrix. The results of this fit are given in the 
last column of Table [U 

We further combine the results of this analysis with 
B — ► J)(*)+/°D^*- ) exclusive branching fractions from 
[H El El El and the BABAR results for B(B — > 
D sJ (2460)-D^) 0, obtaining the following branching 
fractions: 

B(D.j(2460)- -> D* s -n a ) = (56± 13 sta t. ±9 syst .)%, 
S(£> s j(246Q)- -+ Dji) = (16 ± 4 8tat . ± 3 syst .)%, 
B{D- -► 4m-) = (4.62 ± 0.36 stat . ± 0.50 syst .)%. 

In conclusion, we have measured the branching frac- 
tions for the decays B -> dW + ' d|* )_ . These are con- 
sistent with the existing measurements and, in several 
cases, have a significantly smaller uncertainty. The com- 
bination of these results with the existing measurements 
provide the branching fraction for D~ — > 4>tt~ , which 
is also consistent with the most recent measurement 
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TABLE I: Event yields (Nat), efficiencies (e), and branching fractions (B) for pairs of detected decay modes, separately 
and combined. In this combination we use only the results in this paper. B(D7 — > <M~) is a f ree parameter and is also 
reported in the table. The first uncertainty on B is statistical, the second is systematic. The parameter k corresponds to k = 
3.6%/(B(D- 07T-)). 



Decay mode 


-^nicas 


Nat 


e{%) 


B{%) 




Combined B(%) 




, n- n+ 
> u s u 




D+ 
DJ 


86 ±17 
39 ± 9 


3.29 ±0.16 
1.79 ±0.12 


0.90 ±0.18 ±0.14 
(0.74 ±0.17 ±0.13) ■ 


k 


0.64 ±0.13 ±0.10 




, n* _ n+ 




D+ 
D*- 


63 ±19 
30 ± 9 


3.24 ± 0.16 
0.91 ±0.08 


0.67 ±0.20 ±0.11 
(1.15 ±0.33 ±0.26) • 


k 


0.69 ±0.16 ±0.09 


B° 


— > n~ n*+ 




D *+ 
D7 


48 ±13 
68 ± 12 


2.86 ±0.13 
1.63 ±0.10 


0.57 ±0.16 ±0.09 
(1.42 ±0.26 ±0.20) ■ 


k 


0.71 ±0.13 ±0.09 


R° 
13 






D*+ 

or 


129 ± 18 
84 ± 14 


2.68 ± 0.09 
0.86 ± 0.05 


1.65 ±0.23 ±0.19 
(3.38 ± 0.60 ± 0.61) ■ 


k 


1.68 ±0.21 ±0.19 


B~ 






D° 
D7 


214 ±28 
66 ±10 


3.46 ± 0.11 
1.28 ±0.07 


1.33 ±0.18 ±0.32 
(1.11 ±0.17 ±0.17) • 


k 


0.92 ± 0.14 ± 0.18 


B~~ 






D° 

or 


160 ±31 
26 ± 10 


3.71 ± 0.12 
0.64 ± 0.05 


0.93 ±0.18 ±0.19 
(0.87 ±0.33 ±0.16) ■ 


k 


0.77 ±0.15 ±0.13 


Br 






D*° 
D7 


152 ±29 
52 ± 11 


2.69 ±0.10 
1.33 ±0.07 


1.21 ±0.23 ±0.20 
(0.82 ±0.18 ±0.10) ■ 


k 


0.76 ±0.15 ±0.13 


B' 






D*° 
D*~ 


216 ±33 
90 ±15 


2.73 ± 0.07 
0.82 ± 0.04 


1.70 ±0.26 ±0.24 
(2.38 ±0.41 ±0.31) ■ 


k 


1.62 ±0.22 ±0.18 


D7 


— > 07T~ 








4.58 ± 0.48 ± 0.68 


B° 


-> A,j(2460)" 


D+ 


D+ 


27 ± 16 


3.61 ±0.27 


0.26 ±0.15 ±0.07 






B° 


-> D s j(2460)" 


D*+ 


D*+ 


64 ± 15 


2.51 ±0.15 


0.88 ±0.20 ±0.14 






B- 


-> D sJ (2460)- 


-D° 


D° 


75 ±28 


3.78 ± 0.24 


0.43 ±0.16 ±0.13 






B~ 


D sJ (2460) 


-D*° 


D*° 


147 ± 34 


2.81 ±0.14 


1.12 ±0.26 ±0.20 







jlfij and confirms a larger value compared to the pre- 
vious world average 0. We have extracted the absolute 
branching fractions for B -> D(*)+'°D s , / (2460)", thus 
allowing the first measurement of the D s j(2460)~ decay 
rates. Our results show that the D s j (2460) ~ meson de- 
cays via photon or tt° emission to Ds in (72 ± 19)% of 
the cases. 
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